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I. 


INTRODUCTION 


This  report  has  been  prepared  as  an  extension  of  the  work  reported 
earlier,*  to  apply  the  same  techniques  to  a second  set  of  turbulence-related 
problems  in  laser  semi-active  terminal  guidance.  The  problems  of  interest 
here  concern  the  following: 

1)  The  expected  size  of  the  laser  spot  at  the  target. 

2)  Estimation  of  propagation  conditions  from  the  laser  designator 
position. 

3)  Calculation  of  the  power  spectrum  of  the  fluctuations  of  the 
apparent  position  of  the  laser  spot,  as  seen  by  the  laser-spot 
tracker. 

These  matters  shall  be  discussed  in 'the  following  sections.  In  general,  a 
notation  that  is  compatible  with  that  of  the  previous  work  will  be  used  and  an 
awareness  on  the  part  of  the  reader  with  the  notation  presented  there  will  be 
assumed. 


II.  LASER  DESIGNATOR  SPOT  SIZE 

Laser  designator  spot  size  is  determined  by  three  factors.  These 
are  related  to  the  following: 

1)  The  diffraction-limited  beam  spread. 

2)  The  turbulence-induced  beam  spread. 

3)  The  deliberately  induced  transmitter  beam  spread. 

The  first  two  of  these  are  most  conveniently  treated  in  combination,  making 
use  of  the  turbulence-limited,  short-exposure  resolution  theory  as  developed"* 
and  applied.*  The  effective  beam  spread  due  to  these  two  sources  is  shown  to 
be  given  by  the  expression 

♦Fried,  D.L.,  Computer  Simulation  of  Turbulence-Induced  Pointing  Jitter  for 
a Laser  Designator.  Optical  Science  Consultants  Report  No.  TR-203,  March 
1976. 

** Fried,  D.L.  /'Optical  Resolution  Through  a Randomly  Inhomogeneous  Medium 
for  Very  Long  and  Very  Short  Exposures,"  J.  Opt,  Soc.  Am. . Volume  56, 
1372,  1966,  p.  11,  Equation  (18). 

+Frled,  D.  L. , Theoretical  Study  of  Non-Standard  Imaging  Concepts.  Final 
Technical  Report  on  Contract  No.  F30602-74-C-0115,  Rome  Air  Development 
Center  Report  No.  RADC-TR-76-51. 
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The  quantity  rQ  is  defined  by  the  expression 


= { 16.70  X"2  / dsC  2(s/R)5/3l 

' Path  ‘ 
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(4) 


where  the  integration  path  is  from  s = 0 at  the  laser  designator  to  s = R at  the 

2 

designated  target  range.  CN  is  the  refractive-index  structure  constant,  X is 

the  laser  wavelength,  and  D Is  the  diameter  of  the  laser  designator  aperture.* 

80  may  be  taken  to  represent  a full  width  at  half  maximum  due  to  diffraction 

l,  i 

and  turbulence  beam  spreading. 

To  combine  this  spread  with  that  due  to  the  deliberate  defocus  of  the 
designator  unit,  where  this  spread  is  to  a full  width  at  half  maximum  of  6 0 , 

u 

a quadrature  combination  is  utilized.  Thus  the  following  equation  can  be  written 
for  the  laser  spot  angular  spread  due  to  all  three  causes: 


603  = l<60l,2)2+(603)2ll/2 


(5) 


*Op.  Cit. , Computer  Simulation  of  Turbulence-Induced  Pointing  Jitter  for 
a Laser  Designator. 


At  a range  R,  this  corresponds  to  a spot  diameter 


6X3  = R 


(6) 


This  is  the  size  of  the  laser  spot  as  projected  on  the  target. 

III.  ESTIMATION  OF  PROPAGATION  CONDITIONS 

To  form  an  estimate  of  the  propagation  conditions,  it  is  suggested 

that  the  operator  of  the  laser  target  designator  view  the  target  with  high  power 

optics  (binocular-type  optics)  and  try  to  find  an  edge  or  line  on  the  target  from 

which  he  can  determine  the  turbulence-limited  angular  resolution.  It  is 

hypothesized  that  a trained  operator  will  be  able  to  judge  by  comparison  the 

sharpness  of  such  details  with  a set  of  reticle  line-width  markings  what  the 

turbulence-limited  resolution  is,  and  thus  what  the  value  of  r is.  In  practice, 

o 

he  would  simply  be  concerned  with  angular  resolution,  but  it  is  more  convenient 
to  assume  that  he  is  evaluating  an  r related  quantity.  This  quantity  is  denoted 
by  S and  defined  as  ° 


S = / ds  C 2 (R  - s)/R 

Path 


(V 


where  s = 0 is  at  the  observer/laser  designator  end  of  the  propagation  path  and 
s = R is  at  the  target  end. 

The  smaller  S is,  the  better  the  propagation  conditions  are.  It  is 
suggested  that  the  quantity  S be  evaluated  for  all  simulations  and  provided  as 
a printout  associated  with  the  simulation.  It  is  expected  that  to  the  extent  that 
turbulence  conditions  influence  laser  target  designation  performance,  the 
variations  in  $ will  be  correlated  to  the  variations  in  performance. 

It  should  be  recognized  that  this  hypothesized  relationship  is  only  a 
plausible  conjecture  and  presumes  a certain  well-behavedness  of  the  turbulence 
conditions.  For  example,  the  range  weighting  in  evaluating  S is  the  reverse  of 
that  for  the  laser  beam  spread  so  that  the  relationship  would  be  expected  to  hold 

2 

only  if  the  values  of  C^T  are  reasonably  symmetrically  distributed  along  the 

propagation  path.  The  use  of  & as  a measure  of  the  propagation  conditions 
ignores  real  and  pseudo  wind  velocity.  This  is  based  on  the  assumption  that 
normal  winds  will  be  sufficiently  fast  that  the  operator  will  not  be  able  to  track 
out  turbulence-induced  pointing  errors. 


Assuming  that  the  preceding  conjecture  regarding  the  relationship  of  the 
value  of  8 to  the  turbulence  effects  on  designator  performance  is  validated  by 
the  computer  simulations,  it  would  then  be  necessary  to  demonstrate  that  an 
observer  can,  with  as  little  equipment  as  a reticle  in  a pair  of  binoculars,  form 
a reasonably  accurate  estimate  of  8 . This  sort  of  thing  has  been  done  by 
Lincoln  Laboratory  viewing  a specially  prepared  target  board.*  Field  tests 
will  have  to  be  performed  to  see  if  acceptable  results  can  be  obtained  using 
unprepared  targets. 

IV.  ANGLE-OF-ARRIVAL  FLUCTUATION  POWER  SPECTRUM 

In  addition  to  the  turbulence-induced  wander  of  the  laser  spot  on  the 
target,  there  is  a second  source  of  apparent  target  position  jitter  due  to  the 
apparent  variations  of  the  laser  spot  position  associated  with  the  viewing  pro- 
cess. Even  if  the  laser  spot  were  perfectly  stationary,  its  apparent  position 
would  fluctuate  due  to  the  turbulence  in  the  path  from  the  laser  spot  to  the 
tracker  optics.  Associated  with  this  process,  there  is  a viewing  process 
fluctuation  power  spectrum,  Fv(f) , which  should  be  added  to  the  laser  spot 

position  power  spectrum,  (f) , calculated**  to  yield  the  total  power  spectrum, 

FT(f) , associated  with  the  turbulence-induced  jitter  between  the  laser  tracker 

reported  target  position  and  the  actual  target  position. 

FT(f)  = Vf)  +FV(f)  ’ (8) 

It  is  F^,(f)  which  determines  turbulence  contribution  to  the  error  in  a laser 
target  designator/laser  semi-active  seeker  missile  system. 

For  clarity  of  notation,  Dy  shall  be  used  to  denote  the  laser  seeker 

optics  diameter,  and  R^.  to  denote  the  range  from  the  seeker  to  the  target. 

Following  the  same  analysis**  the  nominal  spectrum  for  F^(f)  can  be  written 
as  follows: 


*The  Lincoln  Laboratory  experiment  used  a v-shaped  pattern  of  pairs  of  lights; 
the  observer  noted  how  close  to  the  vertex  of  the  v he  could  resolve  the  pair 
of  lights.  He  used  this  to  infer  rQ  and  could  have  inferred  the  value  of  8 . 

**  Op.  Cit.  Computer  Simulation  of  Turbulence-Induced  Jitter  for  a 
Laser  Designator. 
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(9) 


Nv 

V£>  = l Fv,i(f) 


where 


FVjl(f)  - 1.32xl0-2(x/Dv)2(Dv/r0iVil)5/3F0iVil-1/3f-2/3Gv,f/f0iVil) 


The  function  G ( f/ f . ) is  defined  by  the  expression 

V V)  1 


if  0 < f < 0. 332  f. 


0,V,i 


GV(f/f0ti)  =<  1*12-0.361  (f/f0tV>i)  if  0.332  f0tVf.<fs3.10fo^. 


i£  3*10£o,v,i<f 


(11) 


Here  the  quantities  r.  „ . and  F 


0,V,  i 0,  V,  i 


. are  determined  by  the  local  turbulence 
* th 

and  effective  wind  characteristics  in  the  i segment  of  the  propagation  path 
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from  the  laser  spot  to  the  seeker  optics.  With  C denoting  the  refractive- 

th  N,V,i 

index  structure  constant  in  the  i segment,  V „ . denoting  the  effective 

th  ett*  V* 1 

wind  velocity  in  the  i segment,  z denoting  the  position  on  the  z-axis  of  the 
th  * 1 

midpoint  of  the  i segment,  and  A z^  denoting  the  width  of  the  l segment 
the  following  equation  can  be  written: 


0,V,i 


ll6-,£‘2v,1CN.v„‘V,./V5/V'‘3/5 


(12) 


and 


eff,  V.  i 


O.V.i  -Dv(zv>i/Rv) 


(13) 


.-%*»...■»» 


It  is  to  be  recalled  that  z . is  smallest  near  the  laser  spot  and  largest  near 
the  seeker  optics.  * 1 

The  calculation  of  F^( f)  is  straightforward  and  follows  directly  from  the 

corresponding  sample  calculation  presented.*  In  the  case  under  consideration 
here,  however,  there  is  one  more  effect  to  be  considered  which  makes  it 
necessary  to  distinguish  between  F and  F^.  This  additional  effect  has  to  do 

with  angle-of-arrival  isoplanatism.  Because  of  the  size  of  the  laser  designator 
spot,  fix^,  the  turbulence-induced  angle-of-arrival  has  different  values  for 

light  arriving  from  different  portions  of  the  laser  spot.  This  is  due  to  a lack  of 
angle-of-arrival  isoplanatism. 

To  accommodate  the  lack  of  angle-of-arrival  isoplanatism,  it  is  suggested 
that  the  spot  be  considered^ to  be  subdivided  into  a set  of  N independent  source 
regions  for  each  of  which  F^(f)  applies,  but  for  each  of  which  the  angle-of- 

arrival  fluctuations  are  independent.  Then,  to  a first  approximation,  the 
following  equation  can  be  written: 


Fy(f)  = Fy(f)/N 


It  can  be  seen  that  this  formulation  is  only  approximate.  It  is  expected  that  the 
lower  frequency  angle-of-arrival  fluctuations  be  correlated  over  a larger  por- 
tion of  the  laser  spot,  while  the  higher  temporal  frequency  fluctuations  would 
be  correlated  over  a smaller  portion  of  the  laser  spot.  To  take  this  into  account 
property  would  involve  a very  complex  calculation  to  develop  a value  for  N as  a 
function  of  temporal  frequency  f. 

This  complexity  is  avoided  by  arguing  that  it  will  be  sufficient  if  the  single 
value  of  N is  chosen  so  that  the  resultant  power  spectrum  Fy(f)  will  correspond 

to  the  true  variance  of  the  effective  angle-of-arrival,  i.e. , angle-of-arrival 
averaged  over  the  total  laser  spot.  In  this  way,  there  is  assurance  that  although 
there  may  be  some  distortion  of  the  power  spectrum  of  the  effective  angle-of- 
arrival  for  viewing  the  laser  spot,  the  power  spectrum  will  yield  the  correct 
total  variance  of  the  angle-of-arrival. 

To  determine  N,  an  angle-of-arrival  isoplanatism  angle,  i!>  must  be 
defined.  This  angle  is  defined  as  the  angular  separation  of  two  point  sources 

♦Ibid. 


such  that  the  correlation  of  angle-of-arrival  fluctuations  for  the  two  sources  is 
just  one-half  of  the  variance  of  the  angle-of-arrival  for  either  source.  Then  it 
may  be  considered  that  there  are 


independent  regions  on  the  laser  spot.  The  problem  at  this  point  clearly 
l'educes  to  that  of  determining 

The  subject  of  angle-of-arrival  isoplanatism  has  previously  been  studied.* 
It  has  been  shown  that  the  mean  square  difference  in  angle-of-arrival  from  two 
point  sources  separated  by  an  angle  ^ can  be  written  as** : 


F ( & s/D) 
a 


(16) 


*Op.  Cit. , Theoretical  Study  of  Non-Standard  Imaging  Concepts,  p.  107  et  seq. 

**Ibid. , p.  135,  Equation  (2). 


t 9 
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As  presented,  Equation  ( 1G)  applies  for  an  infinite  plane  wave  source  and 
s = 0 corresponds  to  the  location  of  the  sensor  optics  of  diameter  D.  To  take 

account  of  the  point  source  nature  of  the  matter,  a factor  of  | (R  -s)  'R  | ° ** 

V V 

is  introduced  into  the  integrand*  in  Equation  ( 16).  Moreover,  to  maintain  con- 
sistency with  our  other  formulations,  the  variable  of  integration  is  reversed  so 
that  s = 0 corresponds  to  the  laser  spot  position  and  s R r corresponds  to  the 

sensor  optics  position.  Moreover,  taking  note  of  the  fact  that  the  sensor  optics 
diameter  is  denoted  by  for  this  problem,  the  following  equation  can  be 

written  in  place  of  Equation  ( 16): 


| 


r 


ta*.  

. > ■ — » y v 1 


Dl,a(^  = 


D 


-1/3 


I 

Path 


ds  C 


N 


(s/Ry) 


5/3 


U(Rv-s)/Dy| 


(20) 


The  evaluation  of  F (^)  is  a straightforward,  if  tedious,  numerical  cal- 

O' 

culation  based  on  Equations  (17),  (18),  and  ( 19) . These  calculations  have  been 
carried  out  and  the  numerical  results  listed  in  the  table  were  obtained.  In 
terms  of  a computer  calculation,  it  is  suggested  that  interpolation  on  the  data 
in  the  table,  rather  than  performance  of  the  two-dimensional  integration  of 
Equation  ( 17) , is  the  preferred  approach.  In  either  case,  F(  £)  is  considered 
to  be  a well-defined  function  and  attention  may  be  redirected  to  the  calculation 
of  the  angle-of-arrival  isoplanatism  angle,  d ^ . 


The  mean  square  difference  in  angle-of-arrival  as  a function  of  angular 
separation  is  equal  to  twice  the  angle-of-arrival  variance  minus  twice  the 
angle-of-arrival  covariance,  i.  e. , 


-C,  (d) 

l.o  ' 


(21) 


2 

where  <j  is  the  angle-of-arrival  variance,  and  C,  (d)  is  the  covariance  as 
a 1 , ( v 

a function  of  angular  separation,  d . When  d is  very  large,  the  covariance 
vanishes  and 


lim 

ij—oo 


D , (d) 

1,0 


(22) 


*Ibid. , p.  120,  Equation  (19). 
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Thus,  the  angle-of-arrival  isoplanatism  angle,  «s>  , can  be  defined  for  which  the 
covariance  is  equal  to  one-half  the  variance,  i.  e. , 


C1  'V  = 2 a 2 

l,o  0 2 o' 


by  the  expression 


(J  ) = -D  (a.) 
l,o  0'  2 l,o  ’ 


Making  use  of  tlie  data  in  the  table,  the  following  equation  can  be  written  from 
Equation  (20): 


Dy  (°°)  = 11.  97  D 

l,o  \ 


I dsCN2  (s/V1 

Path  v 


With  the  value  of  D (“)  thus  specified,  it  is  then  a straightforward  matter 

i,Q’ 

to  iterate  on  the  evaluation  of  ( & ) as  defined  by  Equation  ( 20)  to  find  a 

value  of  & for  which  D (i?)  is  one-half  the  value  of  D , (“)  as  calculated 

1,0'  l,o  ' 

from  Equation  (25).  This  value  of  $ corresponds  to 

\\  ith  evaluated,  Equation  ( 15)  may  be  used  to  calculate  N.  The 

evaluation  of  Fy(f)  from  Fv(f)  is  then  a trivial  matter.  The  calculation  of 

F yl  f)  has  been  set  up  in  Equations  ( 9)  to  ( 13) . The  final  evaluation  of  Ft f t) , 

the  total  observed  laser  spot  position  fluctuation,  including  both  transmitter 
and  receiver  jitter,  incorporating  the  calculations*  is  now  a straightforward 
matter. 


*Op.  Cit. , Computer  Simulation  of  Turbulence-Induced  Pointing  Jitter  for  a 
Laser  Target  Designator. 
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TABLE.  CALCULATED  VALUES  OF  F (4) 

a 


The  values  of  £ in  the  table  are  logarithmically  uniformly  spaced,  four  to 

1/4  2/4  *1/4 

decade,  (viz.  1.778  = 10  , 3.162  = 10  , 5.623=  10  ' ). 


f* 


£ 

V<> 

-3 

-6 

1.000  x 10 

1.  069  x 10 

1.778 

3.376 

3.  162 

1.067  x 10~4 

5.623 

3.  376 

-2 

-3 

1.  000  x 10 

1.067  x 10 

1.778 

3.  377 

3. 162 

1.  066  x 10-2 

5.  623 

3.  356 

1.049  x 10_1 

1.778 

3.204 

3.162 

9.205 
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